arXiv:1502.03920vl [astro-ph.HE] 13 Feb 2015 


Slowly rotating superfiuid neutron stars with isospin dependent entrainment in 

a two-fluid model 


Apurba Kheto and Debades Bandyopadhyay 
Astroparticle Physics and Cosmology Division and Centre for Astroparticle Physics, 

Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata-700064, India 

Abstract 

We investigate the slowly rotating general relativistic superfluid neutron stars including the entrain¬ 
ment effect in a two-fluid model, where one fluid represents the superfluid neutrons and the other is the 
charge-neutral fluid called the proton fluid, made of protons and electrons. The equation of state and the 
entrainment effect between the superfluid neutrons and the proton fluid are computed using a relativistic 
mean field (RMF) model where baryon-baryon interaction is mediated by the exchange of a, uj, and p mesons 
and scalar self interactions are also included. The equations governing rotating neutron stars in the slow 
rotation approximation are second order in rotational velocities of neutron and proton fluids. We explore 
the effects of the isospin dependent entrainment and the relative rotation between two fluids on the global 
properties of rotating superfluid neutron stars such as mass, shape, and the mass shedding (Kepler) limit 
within the RMF model with different parameter sets. It is observed that for the global properties of rotating 
superfluid neutron stars in particular, the Kepler limit is modified compared with the case that does not 
include the contribution of p mesons in the entrainment effect. 

PACS numbers: 97.60.Jd, 47.75.-|-f, 95.30.Sf 
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I. INTRODUCTION 


The study of superfluid dynamics in neutron stars has gained momentum recently with the 


observation of fast cooling of the neutron star in Cassiopeia A (Cas A) [l|. It has been inferred 
that the rapid cooling in the neutron star in Cas A might be the outcome of neutron superfluidity 
in its interior^. The glitch phenomenon in neutron stars might also be strong evidence of the 
superfluidity in the crust and core of a neutron star j^b]. There might be an interplay between the 
superfluidity in neutron stars and superfluidity studied in the laboratory. One important aspect 
of the superfluidity is the entrainment effect which was found in a mixture of superfluid ^He and 
'^He in the laboratory Q]. A similar effect might occur in superfluid neutron stars when neutron 
and proton fluids pass through each other. In this case, the two fluids are coupled because the 
momentum of one fluid carries along with it some mass current of the other fluid. This is known 
as the entrainment effect. 

The entrainment effect in superfluid neutron star matter was calculated in relativistic mean field 


(RMF) models 


0 , 0 . 


Comer and Joynt exploited the a-co Walecka model for this purpose. However, 
neutron star matter is highly asymmetric, and the role of symmetry energy is very important in 
determining the equation of state (EoS) and the structures of neutron stars. It is expected that 
the symmetry energy might also influence the entrainment effect. Recently we investigated the 
entrainment effect in the the RMF model including p mesons We showed that the symmetry 


energy signihcantly affected the entrainment effect compared to the case without p mesons [9[. It 
may be worth mentioning here that the dependence of the entrainment effect on the symmetry 


energy was also studied using polytropic equations of state 
liquid theory 12l-ll5l|. 


Q[n 


111 ] as well as with relativistic Fermi 


The role of the entrainment effect in rotating neutron stars was investigated in Newtonian as well 
as general relativistic formulations by different groups [l^, 11, Q . In some of those calculations, 
the dependence of the entrainment effect on the symmetry energy was considered through the 
polytropic EoS [l^, 111. However, so far, there is no calculation of rotating neutron stars based on 
the isospin dependent entrainment effect derived from a realistic EoS. 

In this paper, we are interested in the role of isospin dependent entrainment on slowly rotating 
superfluid neutron stars. Here we adopt the two-fluid formalism for slowly rotating superfluid 
neutron stars as described in Ref, [l^ . The paper is organised in the following way. In Sec. H 
we describe the formalism for calculating the isospin dependent entrainment in a RMF model of 
dense baryonic matter and the application of Hartle’s slow rotation approximation to Einstein’s 
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field equations for superfluid neutron stars. We discuss results in Sec. III. Section IV gives the 
summary and conclusion. 


II. METHODOLOGY 


A. The superfluid formalism 


Here we consider the superfluid formalism developed by various groups 17 h 23I|. The signature 
of the metric used here is the same as in Ref.Q]. The master function (A) in the superfluid 
formalism is a function of three scalars, = —n^n^, = —p^p^ and = —nfj_p^ which are 
constructed from neutron (n^) and proton (p^) number density currents. It may be noted that 
—A(n^,p^,x^) corresponds to the total thermodynamic energ y d ensity when neutron and proton 
fluids are comoving. The stress-energy tensor is written as [Sl. Il6 


+ p^^xv 


( 1 ) 


and the generalized pressure is given by 


= A - nPpp - p^Xp ■ 

The neutron and proton momentum covectors 


Pu = Briu + Apu , 
Xu = An^ Cpu , 


( 2 ) 

( 3 ) 

( 4 ) 


are conjugate to and p^, respectively. It is manifestly evident that neutron or proton momentum 
is a linear combination of both number density currents. The magnitudes of which are chemical 


potentials of neutron and proton fluids, respectively [s]. The charge-neutral proton fluid is com¬ 
posed of protons and electrons, and it was shown that the chemical potential of the proton fluid is 

nn 

the sum of proton and electron chemical potentials [2^, [2^ . The master function is independent 
of the entrainment effect; i.e., = 0 when the coefficient A is zero. One obtains the coefficients 

of Eqs. (I3|) and Q from the master function. 


. dA r, ^dA 


( 5 ) 


The field equations for neutrons and protons involve two conservation equations as well as two 
Euler equations. 
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In the slow rotation approximation, the master function is written in terms of — np, which 
is small with respect to np 


A(n 2 ,p^x 2 ) = -np)* _ 


( 6 ) 


i=0 


Using this form of master function, the coefficients A, etc., that determine the nonrotating 

nn 

background configuration are calculated easily [ 8 |, [ 1 ^ . 

For the slow rotation approximation, we are interested in terms up to second order in the 
rotational velocities of neutrons and protons. This corresponds to the terms proportional to — np 
in the master function. It may be noted that the following combinations appearing in the field 
equations are dependent on 71 when computed on the background [l^ : 


dA dA 971 { d^i\ / 2 

A + n— + np—^ = -71 - n— -> . 7i + n— (x - np) 

on ox^ on \ on ' 

i=2 ^ ' 

,, dA ^ f , d-fA , 2 N*-i 

A + p— + npj^ = -71 - p— - 2 ^ 7 i + (a- 

Qp nr) \ ^ 


dp 


i=2 


dp J 


— npY 


( 7 ) 

( 8 ) 


The calculation of the master function in the RMF model was described in detail in Refs. I 

Q 

Unlike the calculation of Comer and Joynt [ 0 ], the role of symmetry energy on the master function 
and the entrainment effect was considered in Ref. [911 . In the latter case, the relativistic a-to-p 

n Q 

model including scalar meson self-interactions [26|] , was used to derive the master function [9( . The 
Lagrangian density for nucleon-nucleon interaction has the form 

(iltidf" - ms + - gpBlAB ■ 

B=n^p 

+ \ {df^crd^a - mla^) - ^bm {g^af - {gao)'' 
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\pau ■ + Irnlp ■ fA 


( 9 ) 


The Dirac nucleon effective mass m* is defined as m* = m— < g^cr > where the nucleon mass (m) 
is taken as the average of bare neutron (nin) and proton (nip) masses. The frame in which neutrons 
have zero spatial momentum and protons have a wave vector kp = (fco, 0, 0,K) [s] is chosen to solve 


the equations of motion for meson fields in the mean field approximation 


27|. 


The master function, generalized pressure, and chemical potentials of neutron and proton fluids 
in the limit it' —>■ 0 are given by 

.2 _ .2 
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I'-lo ~ 

7 r 2 

U 2 

dA 

dkn 

0 



-kl) + 

\/kl 

+ "i*lo > 

(11) 

Xio = 

7 r 2 

jU2 

rip 

dA 

dkp 

0 

+ kl) + 

c2 

P (k^ 

127 r 2 ^ P 

-kl) + 

yjkj 

+ "i*lo + \/kp+ml , 

(12) 


'^lo ~ ^lo + (a^Io + xio ^p) ) (13) 

where the subscript ”0” stands for quantities calculated in the limit iC —>■ 0, ■, 

4 = [aujlmujf, and = {gp/nipf, and 

m^lo = m^.{kn,kp,0) 


= m — m 


Cg 

*10 2 ,r 2 


kn^kl + m2|p + kp^Jk‘^ + m^lg + ^ m^|gln 


-kn + y/kl + ml\ 


kn + ^kl + 
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. In 
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*10 


K + \ kl + ml\ 


+ bmc^ (m — m^y + cc^ (m — m*)" 


(14) 


T y f^p w Mt^iQ 

It is to be noted here that electrons are treated as noninteracting relativistic particles and are 
included in the calculation of the master ’ function and generalised pressure. The values of the 
various coefficients .A|o, ^lo, C|o) -^olo) l^olo and Cq|o that appear in the field equations are provided 
in the Appendix. 


B. Slowly rotating superfluid neutron stars 


Andersson and Comer [17[ extended Hartle’s slow rotation formalism for the single fluid 


28| 


to the case of the two-fluid model in order to describe superfluid neutron stars. They considered 
that the superfluid neutron and the proton fluid are rotating with different rotational velocities. 
However, they did not include the entrainment effect in their calculation. Here we adopt the 
two-fluid formalism of Andersson and Comer as described by Refs. [16|, llTl] to study stationary, 
axisymmetric, and asymptotically flat configurations. Furthermore we introduce the isospin de¬ 
pendent entrainment in this calculation. In the slow rotation approximation, rotational velocities 
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of neutron and proton (Qp) fluids are considered as small so that inequalities « c and 
^IpR « c are satisfled, where c is the speed of light. The slow rotation acts as the perturbation 
on nonrotating conflgurations. We retain terms up to second order in the angular velocities of 
neutron and proton fluids in fleld e quat ions in the slow rotation approximation. The metric used 
here has the following structure QHQ: 


gp^dx^dx^' = —{N'^ — sin^9K[N'^]'^)dt'^ + Vdr"^ — 2KN'^sin^Odtdcj) + K (d0^ + sin^0d0^) . (15) 


The equations relevant for the metric variables in the two-fluid model and the slow rotation ap¬ 
proximation are same as those of Hartle’s single-fluid model and the metric functions are expanded 


in powers of angular velocities [l6. 17. 


V = (1+ 2u(f,0)) , 


K = r^(l + 2k(f,9)) , 


= uj{r) , 


(16) 


where cu is a first order quantity in angular velocities, and h, v, and k are second order quantities. 
Further h, v, and k are decomposed into i = 0 and £ = 2 terms after expanding those in spherical 
harmonics, 

h = ho{r) + / 12 (f)P 2 (cos 0 ) , 

V = vo{r) + V2{f)P2{cos6) , 


k = k2{f)P2{cos9) , 

where P2(cos0) = (3cos^0 — l)/2. 

Similarly, neutron (n) and proton [p) number densities are expanded as 


(17) 


n = no(f) (1 -I- r]{f, 6)) 
where terms r] and $ are of 0{kl‘^ p), 


p = Po{f){l + <^>{k,9)) , 


(18) 


9 = Vo{r) + r]2{r)P2{cos9) , = <l>o(f) + $ 2 (f)-P 2 (cos 6 ') . (19) 
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A coordinate transformation f —>■ r + C{r,9) is introduced such that A{f{r,9),9) = Ao(r) jl7l |. 
Here the ^ coordinate is also expanded in spherical harmonics as ^ = ^o(^) + C2(^)^’2(cos0). 

With this prescription of the slow rotation approximation for metric functions as well as neu¬ 
tron and proton densities along with the coordinate transformation, the fluid and Einstein field 
equations are reduced to four sets of equations. The hrst set of equations corresponds to nonro¬ 
tating background co nfig urations that are obtained from the solutions of two background metric 
components A and v |l6l. Il7|. Those are given in terms of coefficients of fluid equations, 

^olo^'o +-®olo ^0 + 2 ~ 0 > C'o IqPo + ^olo ^0 + 2 aIo ’ (20) 

where prime denotes differentiation with respect to f and Aq|q, HqI^, and CqIq coefficients are 


obtained from the master function and are taken from Ref.[9(]. The regularity condition demands 


• 0 - 


that A, A', v', Uq and Pq vanish at the origin. The total mass of this configuration is 

M = —dvr / Ao(f)f^df . 

Jo 


( 21 ) 


Next, the frame dragging uj{r), which is hrst order in angular velocities of neutron and proton 


huids, is obtained from the following equation 

, dL 


3, Cl 



^4^-{X+u)l2< 


dr 


- 167re(^-")/2 _ Ao) L„ = 167re^^-^^/\opo - Hp) 


( 22 ) 


This equation has the same structure as that of the single huid except for the nonzero term on the 
right-hand side 28|. Here we dehne = u — and Lp = uj — Hp, which represent the rotational 


frequencies as measured by a distant observer. The boundary condition implies that the interior 
solution of a;(r) matches with the vacuum solution 

2 J 


Ln{R) — ~^n + 


R3 


(23) 


where J is the total angular momentum of the system. The derivative of the solution is also 
continuous at the surface [3]- 

The neutron and proton angular momenta, and Jp, respectively, are given by [1 


Jn — 


Stt 


cR 




jJjQTlQLn T -AqTIqPq (Hfi Hp) 


(24) 


and 


J — 

Jp- 3 




dpp4g(A-i.)/2 


XoPoLp + AqUoPo (flp - fin) 


(25) 
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The total angular momentum J is equal to Jn + Jp- 

The last two sets of equations are 0{Q‘^p) equations. One can obtain ,^ 0 ) ho, and vq 

from i = 0 second-order equations, on the other hand, ^ 2 ; ^ 2 , ^ 2 ) ^ 2 ; V 2 , and k 2 follow from 1 = 2 
second-order equations. A detailed discussion of ^ = 0 and 1 = 2 second order equations and 
numerical techniques to solve those equations can be found in Refs.[l^. 17, 28|. After obtaining a 
complete solution in the slow rotation approximation, one can calculate the quadruple moment of 
the configuration and the rotationally induced change of mass as described in Refs.[l6. 17]. Fur¬ 
thermore, the Kepler frequency of the slowly rotating superfluid neutron star is obtained following 


the prescription of Andersson and Comer 


17 |, 


V R3 W> \B?\2M 4R4M2 4 R \ ^ ’ 


where scaling of J = JVlp, 6M = 6M^lp, ^0 = and ^2 = C 2 ^p with Qp is made and 

3(R3_2M3) / 2M\ 3R^-3R^M-2R^M‘^-8RM^+ 6M^ 

" “ 4M3 2RM2(R- 2M) ■ 


(27) 


It is to be noted that the expression for the Kepler frequency in Eq. (12611 differs from that of 
Eq. (77) of Ref.[l^. This difference originates from the factor at the beginning of the third term 
within the second bracket and the term involving ^0 and ^2 in both equations. We discuss this 
issue further in the next section. 

It is worth mentioning here that the model based on the slow rotation approximation is ap- 
pUcable f„. the taeteet obsetved puW as noted by otheta QQ. Howevet, thia apptoxintation 
breaks down near the Kepler limit |10l |. 


III. RESULTS AND DISCUSSION 


Now we discuss the results of slowly rotating superfluid neutron stars. Nonrotating background 
configurations are obtained by solving Eq. (120^ . In this context, we exploit the RMF EoS which 
includes the isospin dependent entrainment effect a. we ttae OL and NL3 pata^etet aeta in thia 
calculation, both of which are listed in Table I. Central neutron number density is an essential 
input for the calculation of the background configurations. The proton number density in the 
background model is no longer a free parameter because the chemical equilibrium is imposed at 
the centre of the star, i.e., /rjo = x|o [3- The chemical equilibrium is established when both 
fluids are corotating. However, the chemical equilibrium does not hold good for different rotation 
rates of neutron and proton fluids {l^ . [^. Masses and radii corresponding to two nonrotating 

















configurations are also recorded in Table I. The chosen background configurations are just below 
their maximum masses Q]. Furthermore, we consider ?7o(0) = 0 and, consequently, ‘ho(O) = 0 in all 
cases. 

As soon as we know the background configuration, we can calculate the frame-dragging fre¬ 
quency from Ea. (|2‘2p . As we are dealing with the two-fluid system, the central value of and 
relative rotation rate are needed to solve Ea. (l22p |l^ . A rescaled Ea. (l22l) with the definition 

of Ln{r) = Lnl^p is solved to determine the frame-dragging frequency for different values of Qp 
using a hxed relative rotation rate. The boundary condition of the problem demands that the 
interior solution matches with the known vacuum solution given by Eg. (1231) . The frame-dragging 
frequency, is plotted as a function of radial distance {r/R) in Fig. 1 for three different relative 
rotation rates. The left panel denotes the GL parameter set and the right panel represents the NL3 
set. The frame-dragging frequency decreases monotonically from the centre to the surface of the 
star for three relative rotation rates in both panels. This feature of the frame-dragging frequency 
is quite similar to the standard single-fluid result Further it is noted that the frame-dragging 
frequency is always higher for larger values of relative rotation rate. 

Now we discuss numerical solutions of different metric functions of the superfluid neutron star 
in the slow rotation approximation . First we solve the £ = 0 equations and determine ‘I’O) 


ho and vq following the procedure laid down by Andersson and Comer [l^. Metric functions ho 
and Vq match with the vacuum solutions at the surface. The metric function uo(r) as a function 
of radial distance is displayed in Fig. 2 for three different relative rotation rates. The left panel 
shows the results of the GL set and the right panel corresponds to those of the NL3 set. It is 
noted that the metric function vq increases monotonically to the surface and matches smoothly 
with the exterior solution. For the NL3 set, the value of this metric function at the surface is 
always higher than that of the GL set. A new metric function mo is defined in terms of vq and 
A as mo = wq/ exp(A). The radial prohle of mo, which merges with the exterior solution at the 
surface, is shown in Fig. 3 for different relative rotation rates. 

_We solve the 1 = 2 equations in a similar way to that used for solutions oi i = Q equations 
A new variable. A: = A:2 -|- /12, is introduced to solve two coupled first-order equations in /12 


and k 2 


28|. This leads to two coupled differential equations in k and / 12 , which are solved using 


the method described by Hartle [281 ]. In Fig. 4 , the metric functions /iq and /12 are plotted as 
a function of radial distance for different relative rotation rates. The results of the GL and NL3 
sets are shown in the left and right panels, respectively. In both panels, the lower three curves 
denote the metric function /iq and the upper three curves imply the metric function h 2 - Figure 5 
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shows the radial profiles of Co{^) (upper curves) and (lower curves), and in Fig. 6, we have 

k = k 2 + h 2 versus r for the GL (left panel) and NL3 (right panel) sets and three different relative 
rotation rates. In Fig. 5, the magnitude of '^ 2 (^) at the surface in the right panel is quite large with 
respect to that of the left panel when the relative rotation rate is larger than 1 and this function 
is directly related to the deformation of the star due to rotation. 

Figure 7 exhibits the variation of rotationally induced corrections to the neutron number density 
no??o (three upper curves) and nQr ]2 (three lower curves) with radius. Similarly, Fig. 8 represents 
the variation of rotationally induced corrections to the proton number density po^o (three upper 
curves) and po ^2 (three lower curves) with radius. In both cases, the left panel denotes the 
results of the GL set and the right panel corresponds to those of the NL3 set. We explore the 
role of symmetry energy on the rotationally induced corrections to the proton number density by 
comparing two cases with (left panel) and without (right panel) p mesons for the GL set in Fig. 
9. For the case without p mesons, we consider a nonrotating configuration that is just below the 
maximum mass neutron star. The mass and radius of this neutron star is 2.33 Mq and 10.96 
km. It is noted that the corrections to the proton number density are significantly modified in the 
presence of p mesons. 

The deformation of a rotating star is obtained in terms of the ratio of the polar and equatorial 
radii. For the slowly rotating star, this is given by ~ 1 + The ratio of polar to equatorial 

radii as a function of relative rotation rate is plotted in Fig. 10 for the GL (solid line) and NL3 
(dashed line) sets. We consider the proton rotation rate to be equal to that of the fastest rotating 


pulsar having spin frequency 716 Hz 


3ll |. The nonrotating situation is achieved when the relative 


rotation rate approaches zero. Furthermore we find that the rotationally induced deformation of 
the star is larger for the NL3 case than the GL case. This deformation increases with increasing 
relative rotation rate. 

As neutron and proton fluids may rotate at different rates, one of them extends beyond the 
other at the equator. The Kepler limit is obtained from the rotation rate of the outer fluid. 
To determine the mass-shedding (Kepler) limit we have to solve the quadratic equation (|26p for 
kip. When > kip the Kepler frequency is determined by the neutrons; for kip > the 
Kepler frequency is determined by the protons. We calculate the Kepler limit in the RMF model 
including p mesons using the GL and NL3 parameter sets for the background configurations of 
Table 1. The mass-shedding (Kepler) limit klx as a function of relative rotation rate is plotted in 
Fig. 11 for the GL set (left panel) and the NL3 set (right panel). We use the radial profiles of 
the entrainment effect in this calculation of Kepler frequency. The results are qualitatively similar 
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to the previous investigation by Prix and collaborators [10] though the authors in the that case 
used some constant values of entrainment. However, our results are different from those of Comer 
[l^ . For > Qp, the Kepler frequency (solid square) approaches a constant value with increasing 
When < 1 , the Kepler frequency (solid circle) monotonically increases with decreasing 

relative rotation rate, as evident from Fig. 11, whereas the opposite scenario was found in the 
work of Comer [l^. On the other hand, Prix et al. [l^ found that the Kepler limit increased 
monotonically as the relative rotation rate decreased. This is quite similar to our results. The 
difference between our results and those of Comer [l^ may be due to different expressions for VLk 
that we have discussed in connection with Ea. (l26p in Sec. II. Furthermore, Comer [l^ calculated 
the entrainment using the the equation of state obtained in the relativistic a-uj model. Without 
p mesons, the effects of symmetry energy on the entrainment was absent. On the other hand, 
we exploit an isospin dependent entrainment effect calculated in the a-oj-p RMF model for the 
determination of the Kepler limit Qj. We compare the Kepler limit calculated in the RMF model 
with and without p mesons for the CL set in Fig. 12. In both cases, we consider nonrotating 
conhgurations that are just below their maximum masses, as noted in Table I and discussed in 
connection with Fig. 7. The solid line denotes the calculation without p mesons and the dashed line 
represents the case with p mesons. Furthermore, solid squares and circles correspond to allowed 
rotation rates of neutron and proton fluids, respectively. It is noted that the two results differ, as 
is evident from the highlighted part of Fig. 12. 


IV. SUMMARY AND CONCLUSIONS 


We have studied the role of the isospin dependent entrainment and the relative rotation rates 
of neutron and proton fluids on the global properties of slowly rotating superfluid neutron stars 
such as the structures and the Kepler limit in the two-fluid formalism. The two-fluid formalism 
of Andersson and Comer [l^ is adopted in our work. The effects of symmetry energy on the 
EoS and entrainment are studied using the a-uj-p RMF model. The symmetry energy significantly 
influences the rotationally induced corrections to the proton number density. It is found that the 
Kepler limit obtained with the isospin dependent entrainment effect is lower than that of the case 
when the isospin term is neglected in the entrainment effect. The behaviour of the Kepler limit as 
a function of the relative rotation rate in our case is qualitatively similar to the results of Prix et 
al. jl^ obtained using the polytropic EoS. The calculation of slowly rotating superfluid neutron 
stars including the isospin dependent entrainment effect in a realistic EoS is the hrst of its kind. 
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V. APPENDIX 


The values of some useful matter coefficients (see y^that are the inputs of field equations are 
the following 
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TABLE I: Nucleon-meson coupling constants corresponding to the GL and NL3 sets are taken from Refs. 


27 


29|. The coupling constants are obtained by reproducing the saturation properties of symmetric nuclear 


matter as detailed in Ref.Q . All the coupling constants are in fm^, except b and c which are dimensionless. 


The background nonrotating configurations as computed by Kheto and Bandyopadhyay [9| are used here . 
The central neutron wave number fc„(0) is given by fm“^. The mass(M) and radius (i?) are in units of Mq 
and km, respectively. 



cl cl cl b c iy{0) 

fcn(O) 

Xp{0) 

M R 

?7o(0) 

GL 

12.684 7.148 4.410 0.005610 -0.006986 -2.38799 

2.71 

0.24 

2.37 11.09 

0.0 

NL3 

15.739 10.530 5.324 0.002055 -0.002650 -2.33319 

2.40 

0.23 

2.82 13.17 

0.0 
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FIG. 1. The frame-dragging frequency uj{r) is plotted as a function of radial distance {r/R) using the GL 
parameter set (left panel) and the NL3 parameter set (right panel) for three different relative rotation rates 
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FIG. 2. The metric function vo{r) is plotted as a function of radial distance {r/R) using the GL parameter 
set (left panel) and the NL3 parameter set (right panel) for three different relative rotation rates 
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FIG. 3. The metric function mo(r) = rno(r)/exp(A(r)) is plotted as a function of radial distance {r/R) 
using the GL parameter set (left panel) and the NL3 parameter set (right panel) for three different relative 
rotation rates fln/flp. 
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FIG. 4. The metric functions ho(r) (three lower curves) and h 2 {r) (three upper curves) are plotted as a 
function of radial distance {r/R) for three different relative rotation rates iln/^p- The left panel shows the 
results of the GL parameter set and the right panel demonstrates those of the NL3 parameter set. 
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FIG. 5. The radial displacement ^o{r) (three upper curves) and ^ 2 {r) (three lower curves) are plotted as a 
function of radial distance {r/R) for three different relative rotation rates iln/^p- The left panel shows the 
results of the GL parameter set and the right panel represents those of the NL3 parameter set. 
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FIG. 6. The quantity k is plotted as a function of radial distance (r/R) using the GL parameter set (left 
panel) and NL3 parameter set (right panel) for three different rotation rates Q,n/VLp. 
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FIG. 7. The rotationally induced corrections to the neutron number density n.o(^)??o(’') (three upper curves) 
and no(r)? 72 (r') (three lower curves) are plotted as a function of radial distance (r/R) for three different 
relative rotation rates Results of the GL and NL3 parameter sets are shown in the left and right 

panels, respectively. 
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FIG. 8. The rotationally induced corrections to the proton number density (three upper curves at 

rjR = 1) and po{r)(j) 2 {r) (three lower curves at r/R = 1) are plotted as a function of radial distance (r/R) 
using the GL parameter set (left panel) and NL3 parameter set (right panel) for three different relative 
rotation rates 
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FIG. 9. The rotationally induced corrections to the proton number density Poif)<Po(j') (three upper curves 
at r/R = 1) and po{r)(j) 2 {r) (three lower curves at rjR = 1) are plotted as a function of radial distance 
(r/R) for the GL parameter set with (left panel) and without (right panel) p mesons for three different 
relative rotation rates 
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FIG. 10. The ratio of the polar to equatorial radii (Rp/Re) is shown as a function of relative rotation rate 
^In/^p corresponding to neutron stars of masses 2.37 Mq with the GL set (solid line) and 2.82 Mq with 
the NL3 set (dashed line), respectively, considering that Vp = Q.p/2'K is equal to that of the fastest rotating 


pulsar having spin frequency 716 Hz 


3ii. 
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FIG. 11. The mass-shedding (Kepler) limit is shown as a function of relative rotation rate for the 

GL set (left panel) and the NL3 set (right panel). The solid squares (green) show the allowed rotation rate 
of the neutron fluid{iln) and the solid circles (red) show the allowed rotation rate of the proton//u*(i(np). 
The Kepler frequency is the largest of the two. 
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FIG. 12. (Color online) The mass-shedding (Kepler) limit is shown as a function of relative rotation rate 
with (dashed line) and without p mesons (solid line) for the GL set. The solid squares show the 
allowed rotation rate of the neutron fluid (fin) and the solid circles show the allowed rotation rate of the 
proton fluid (flp). 
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